Astronomy &: Astrophysics manuscript no. Mel0653 


February 2, 2008 


(DOI: will be inserted by hand later) 





Relativistic Parker winds with variable effective polytropic index 



Z. Meliani^'2, C. Sauty^'^, K. Tsinganos^, and N. Vlahakis^ 



,1,2 



o 
o 

(N 



> 
o 
o 

o 
o 

:^ 

6 



^ Observatoire de Paris, L.U.Th., F-92190 Meudon, France 

^ Universite de Paris 7, APC, 2 place Jussieu, F- 75005 Paris, France 

lASA and Section of Astrophysics, Astronomy & Mechanics, Department of Physics, University of Athens, 

Panepistimiopolis GR-157 84, Zografos, Greece 

Received ... / accepted 23/06/2004 

Abstract. Spherically symmetric hydrodynamical outflows accelerated thermally in the vicinity of a compact 
object are studied by generalizing an equation of state with a variable effective polytropic index, appropriate to 
describe relativistic temperatures close to the central object and nonrelativistic ones further away. Relativistic 
effects introduced by the Schwarzschild metric and the presence of relativistic temperatures in the corona are 
compared with previous results for a constant effective polytropic index and also with results of the classical wind 
theory. By a parametric study of the polytropic index and the location of the sonic transition it is found that space 
time curvature and relativistic temperatures tend to increase the efficiency of thermal driving in accelerating the 
outflow. Thus conversely to the classical Parker wind, the outflow is accelerated even for polytropic indices higher 
than 3/2. The results of this simple but fully relativistic extension of the polytropic equation of state may be 
useful in simulations of outflows from hot coronae in black hole magnetospheres. 
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1. Introduction 

Relativistic outflows are commonly inferred from ob- 
servations of coUimated winds (jets) in Galactic x-ray 
binaries and supermassive black holes in active galac- 
tic nuclei and quasars (Birctta ct al. 1999, Livio 2002). 
Also, observations infer coronae with rather high tem- 
peratures in microquasars (e.g., ICorbel et al. 2003|l and 
AGN (e.g., [Rozanska fc Czerny 2000] and refs. therein). 
Although there are still some ambiguities in the interpre- 
tation of these observations, usually temperatures up to 
lO^K for electrons/positrons and lO^^K for protons are 
usually inferred. Nevertheless, even if such temperatures 
were to correspond to closed field line regions, by using the 
analogy with the solar corona, one may extrapolate that 
such high temperatures could also exist in the open field 
line regions where the outflow is accelerated. The heating 
mechanism to obtain these high temperatures can be of 
magnetic origin | |Heinz fc Begelman 2000| ). Alternatively, 
accreted material falling onto the central black hole may 
be decelerated via a shock, feeding the magnetosphere 
with a hot plasma. Discussions on the possibility of such 
shock formation in an accretion disk that leads to a 
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heated corona have been already given in the literature in 
the so called CENBOL m odel ( Chakrabarti ct al. 19961 
[Chattopadhyay et al. 2004| Das, 2000, and references 
therein). These papers have also discussed how radiation 
and coupling with the photon distribution has to be taken 
into account properly. 

The inevitable result of such hot atmospheres is that 
they expand supersonically and at large distances thermal 
energy is converted to bulk flow kinetic energy. For exam- 
ple, Ferrari et al. (1985) Das (2000) and Chattopadhyay 
et al. (2004) have suggested that outflows can be ther- 
mally accelerated to relativistic speeds. However, such 
treatments assumed the classical polytropic equation of 
state which prevents from studying consistently the rela- 
tivistic temperatures in the corona and the nonrelativistic 
ones farther away. The first effort to use an equation of 
state appropriate for outfiows containing both ultrarela- 
tivistic and classical temperatures has been used only for 
studying spherical accretion flows and only in the adia- 
batic case ( Mathews 1971, Blumenthal fc Mathews 1976| . 
Moreover, Blumenthal & Mathews (1976) discussed the 
topology of the Mach number variation in the solution of 
an adiabatic wind. 

However, as it is well known IjParker 1960|l the topol- 
ogy of the Mach number does not give us information on 
whether the flow itself is accelerated or not. 

In this paper we extend the equation of state used in 
Blumenthal & Mathews (1976) for nonadiabatic flows and 
our goal is to investigate the efficiency of thermal driving 
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under extreme relativistic conditions. For example, in the 
classical polytropic solar wind theory, the polytropic index 
needs to be smaller than 3/2 in order to obtain an acceler- 
ated wind solution IjParker 1960|l . This means that a min- 
imal extension of the corona is required. Here we explore 
how this value is changed under relativistic conditions. In 
order to simplify the study, we shall focus our attention 
on spherical, steady and radial hydrodynamic outflows at 
large distances from a stationary compact object. Also for 
simplicity we shall use Schwarzschild's metric. 

Furthermore, we will not discuss how the pressure and 
the heating may include coupling with the radiation field, 
especially in the case of pair production, as our goal is to 
keep the discussion at a basic level and this has already 
been done elsewhere (e.g., [Chattopadhyay et al. 2004| ). 

In the next section we outline the governing equations 
with particular emphasis to an equation of state appro- 
priate to relativistic temperatures close to the base of the 
outflow and classical ones further away. In Sec. 3 we de- 
rive for a given asymptotic value of the polytropic index 
the range of allowed locations of the sonic surface and also 
the limits of the asymptotic speed. In Sec. 4 via a para- 
metric study we present a comparison of our model with 
the nonrelativistic Parker wind and relativistic winds with 
constant polytropic index. An astrophysical application is 
outlined in Sec. 5 and the results are summarized in Sec. 6. 

2. Basic Equations 

The flow of a relativistic fluid is governed by conservation 
of the number of particles and energy-momentum, 

(nu");, = 0, (1) 
(T°'');b = with T"'' = (?iu;/c2)u"u'' + P5°^ (2) 

where n is the proper number density of particles 
(electrons/protons, or, electrons/positrons) in the co- 
moving frame of the fluid and u"^ = (7c; 71;) is the 
fluid four-velocity, with v the spatial velocity and 7 = 
(l — ti^/c^) the Lorentz factor. 

The proper enthalpy per particle, w, is the sum of the 
proper internal energy per particle (including the proper 
mass), e, and the proper pressure P divided by n. 



w 



P 

e+ — . 

n 



(3) 



The space-time around the central compact object is 
described by using Schwarzschild's metric wherein in order 
to treat analytically the fluid equations the simplest way 
is to express all quantities in a 3-1-1 split space-time. 



where 



r^d9^+r^ sin(?2d^2 



a = W 1 



1 



re 



(4) 



(5) 



is the lapse or redshift factor induced by gravity at a dis- 
tance r from the central object (e.g., a black hole) of mass 



A^^, written also in terms of the Schwarzschild or gravi- 
tational radius tq. 



(6) 



Assuming spherical symmetry such that the velocity 
is purely radial, v = vf, Euler's equation reduces to the 
following single differential equation. 



1 2 ^ dv^ 

2 c"^ dr 



, d In a dP 



dr dr 



^ = 0. 



(7) 



By combining Eqs. Q and lj3Jl, we get that the varia- 
tions of the specific enthalpy and internal energy are 



1 P 
dw = —dP , de — -^dn . 



(8) 



Similarly, the conservation of the number of particles 
for a spherical wind, Eq. (Q, takes the following simpler 
form, 



dr 



{a'fnvr'^) — . 



(9) 



In the following, we shall adopt an equation of state 
appropriate to describe relativistic temperatures close to 
the central object and nonrelativistic ones further away 
and also derive its ultrarelativistic and classical limits. 

2.1. Closure of the fluid equations 

To close the system of equations, (|7I9|I we should specify 
an equation of state that relates the entropy per particle 
s — s{w, P) — s(e, P) to the pressure, P, and the enthalpy 
per particle, w, or the internal energy per particle, e. In 
the literature fe.g.. IMichel 19721 ICh akrabart fet al. 19961 
pDas 1999) various authors have adopted a classical adi- 
abatic or polytropic equation to replace the equation of 
state 



P (X n 



r 



(10) 



where F is a constant. However, in this form F should 
attain different values whether the particles are in a rel- 
ativistic or a classical thermodynamical regime. It is well 
known that for a monoatomic gas, in adiabatic flow the in- 
dex is 5/3 in the classical regime with nonrelativistic tem- 
peratures (fceP <C TOC^), and 4/3 in the ultra-relativistic 
one (^bT ^ mc^). (Here fee is the Boltzmann constant 
and mc^ is the rest mass energy per particle.) Thus we 
can define an effective polytropic index 



- off 



rflnP 
dlnn 



(11) 



which should vary with temperature. In other words, the 
previous equation of state with a constant index cannot 
account for media with a transition from ultrarelativistic 
to nonrelativistic temperatures, as inferred from observa- 
tions of coronae in AGNs. 
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The scalar isotropic pressure of a single perfect fluid is 
given by IjMathews 19711 |Synge 1957| ) 



pi ( e eo 

P = o"fo 

3 \ eo e 



(12) 



where eo = mc^. Its validity for coUisional and collision- 
less fluids is discussed in detail in Blumenthal & Mathews 
(1976) and we refer the reader to section III of their paper. 

By integrating Eqs. ||SJ) and H12|l wc obtain a gener- 
alized equation of state for an adiabatic (ideal) fluid (see 
IBlumenthal fc Mathews l976l Eq. 3.3), 

- eg = nn^ = Kn^"""^ , (13) 

which represents the equation of state for an adiabatic 
flow. Here Fad = 5/3 is the real adiabatic index related to 
the number of degrees of freedom of the particles. 

For studying stellar interiors, a classical trick to study 
different equations of state and their hardness is to use the 
same equation, called polytropic equation of state, with a 
different value of the adiabatic index, assuming the num- 
ber of degrees of freedom of the particles has changed. 
In this case Fad ^ 5/3, but the system still evolves adi- 
abatically. E. Parker to study the solar wind proposed 
IjFarker 1960|l to follow the same path using Eq. IIUI as- 
suming F being any value between 1 and the adiabatic 
value 5/3 to mimic the coronal heating in a simple but hid- 
den way. As explained in Sauty et al. (1999) this is equiv- 
alent to use a generalized enthalpy and internal energy 
which includes the heating. 

Following Parker's initial approach of the polytropic 
solar wind, we assume also here that the enthalpy and the 
particle density are related by Eq. ||SJ| but relating e to 
the particle density via 



e — en = Kn 



,r-i 



(14) 



with again a value of the constant polytropic index F 
which can be anything between 1 (isothermal) and 5/3 
(adiabatic value), or higher in case of energy losses. The 
basic idea is to generalize to heated relativistic winds the 
approach of Blumenthal & Mathews (1976) for relativis- 
tic accretion flows like Parker did with Bondi's accretion 
models. Of course at this stage w (resp. e) is a general- 
ized specific enthalpy (resp. generalized specific internal 
energy) which includes the energy brought via heating 
from the external medium. We recall in Appendix 1X1 how 
to calculate the extra heating implicitly in the case of a 
monoatomic ionized plasma. 

By combining Eq. (|14|l and JS)) we get a generalized 
relation between pressure and density which includes the 
proper mass energy per particle, cq — mc^ , 



P = 



F- 1 



V4 + ' 



(15) 



for F ^ 1. Note that Eq. H15|l is a generalization of the 
usual form of the so-called equation of state with the in- 



clusion of relativistic thermal effects. Thus the effective 
polytropic index (Eq. [TT|l takes the form 



Fcff = F - 



F- 1 



(16) 



2 1 + p' 

where we have introduced a new dimensionlcss function p 



,r-i 



P - 



(17) 



All thermal effects appear in this single function p, hence 
decreasing the number of free parameters in the model. 

Finally, the local sound speed — (3sC can be written 
as a function of p. 



F - 1 



w dn 



1 - 



2 + (3~F)p 
(l+p)[2+(F+l)p] 



(18) 



The logarithmic derivative of the sound speed with density 
is 



dln/32 



(F + 1)(3 - F)p2 + 4(F + l)p + 4F 



dlnp (p+l)[2F+(F + l)p][2+(F + l)p] ' 



(19) 



and, for F < 3, it is positive. Thus, with F < 3 it follows 
that 



V - 1 



(20) 



with the equality holding in the high temperature limit 
(p — > 00). Note that for F < 5/3, we get the usual con- 
dition Pi < 1/3 (jMichel 1972|l . 

2.2. Nonrelativistic and ultrarelativistic limits 

In the limit of a nonrelativistic fiuid, eth ^ eoj we may cal- 
culate the thermal energy per particle or random thermal 
energy IjMathews 1971|l . 



eth 



= e - eo 
K „r-i 



eo-eo 
if eth < eo • 



(21) 



In such a case, p is twice the ratio of the thermal energy 
per particle and the mass energy 



p« 2 — 

eo 



(22) 



For classical temperatures, wherein the thermal energy is 
negligible compared to the mass energy, eth ^ eo, we have 
that p ^ 1 and 



P ~ (F - 1) neth =^ (F-1)— n^, 

zeo 



(23) 



such that the effective polytropic index is exactly Fctr — F. 
In general, as temperature decreases with distance, F al- 
ways represents the asymptotic value of the effective poly- 
tropic index. This allows us to have a simple albeit arti- 
ficial way to model relativistic coronae as in the classical 
Parker wind. 
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On the other hand, in the ultrarelativistic domain 
where the thermal energy is much larger than the mass 
energy (eth ^ cq 4^ 1) we have 

-j^ ^ r 1 

Pc^—^nin^. (24) 

The effective polytropic index becomes now Toff = (F + 
l)/2. 

The effective polytropic index and corresponding tem- 
perature for an adiabatic flow with relativistic tempera- 
tures at the base and classical behavior further away are 
plotted in Fig. These solutions are discussed in more 
detail later on. We simply note here that for an adia- 
batic flow, the polytropic index increases from 4/3 at the 
base where the thermal energy equals or exceeds the mass 
energy (i.e. T > for protons, or, T > lO^K for 

positrons) to 5/3 asymptotically where the temperature 
is much lower (i.e., T <C 10^^ K for protons or T ^ 10^ K 
for positrons), as expected. To emphasize this property of 
the equation of state, we have plotted in Fig.QJa case with 
ultrarelativistic temperatures at the base. Note that the 
coronal temperature close to the compact object are very 
sensitive to the value of the parameter /u, as discussed in 
the following. 




' 2 4 6 2 4 6 

1 10 10 10 1 10 10 10 

r/r r/r^ 

G G 



3. Equations in dimensionless form and 
parameters 

3.1. Momentum and continuity equations 

In the following all quantities are defined in dimensionless 
form. First, velocities are normalized to the speed of light 
P = v/c such that /3s = Vs/c is the (variable) dimensionless 
sound speed. Second, distances could have been normal- 
ized at the gravitational radius tq. However, we found 
easier to normalize all quantities at the sonic surface r^, 
such that we can define a dimensionless radius. 



R = r/r^, . 

A crucial parameter of our model is 



re 



(25) 



(26) 



which is the square of the escape speed at the sonic ra- 
dius normalized to the speed of light. It is also the ratio 
of the gravitational radius to the sonic distance. Thus /j, 
measures the strength of the gravitational field and the 
distance between the central black hole and the sonic sur- 
face.^ 

The Euler equation and the continuity equation can 
be written in dimensionless forms, 



-/3g d\n(3^ 
1-/32 dR 



R 



R R- fi 



(27) 



^ This parameter /i is similar to the parameter m used by 
Daigne & Drenkhahn (2002) except that we use the sonic sur- 
face while they use the Alfven radius for normalization. 



Fig. 1. (a) Plot of the effective polytropic index Feff of 
an adiabatic outflow showing the transition of Fofi from 
the ultrarelativistic value of 4/3 to the classical one 5/3. 
The corresponding temperature profile is plotted in (b) 
both, for electron/proton [k-gT = P/{2n)] and pair plas- 
mas [k^T = P/n] and has ultrarelativistic temperatures 
at the base which decrease to classical ones further out. 



/32 -/32 d\np 



/32 



dR 



(F-1) 



2 
R 



1-/32 1 



/32 2RR- p 



(28) 



In a nonrelativistic wind, the radius of the star is in- 
dependently given from its mass. Conversely if the central 
object is a black hole, the gravitational radius provides a 
natural length scale. However, the wind cannot start and 
accelerate at this particular distance. In fact, it has to start 
at a radius Tq located obviously above the Schwarzschild 
radius in the corona. In other words the wind solution 
should start at a dimensionless radius Ro > p which will 
be an output of the integration starting at the sonic sur- 
face and integrating up wind. 

3.2. Integration of equations and range of values of fi 

This model has two free parameters that affect the solu- 
tion of the differential equations lt?7|l . and 1(^ .2 The first 
parameter, F, is the value of the effective polytropic index 



^ The other parameters are the mass of the central object 
M.*, the mass-loss rate m, and the constant /t. 
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at infinity. Tfie second parameter, fi, is the ratio of the 
Schwarzschild's radius to the radius at the sonic surface. 

Eqs. H27|) and (|28|l have a singularity corresponding to 
the sonic surface, which we may call hydrodynamic hori- 
zon by analogy with the black hole horizon. At the sonic 
surface 

the right hand sides of the two equations vanish be- 
cause the left hand sides vanish. 

For a detailed discussion of the nature and posi- 
tion of the critical surfaces in astrophysical winds, see 



Tsinganos et al. (1996) 



This gives the following criticality condition. 



(29) 



4- 

Note that the value of p at the sonic surface, pi,, follows 
from Eg. lfTH)) . 

Since /3g^ < 1, the above criticality condition requires 
that we should always have /i < 1. Furthermore, by com- 
bining Eqs. (|20() and H29() we find an upper value for p, 



pL < /In 



4(r-i) 
3r- 1 



(30) 



Once we get the value of /3s* (<^ Ps*) from the sonic 
criticality condition, we deduce the values of the slope of 
the velocity at the critical surface using De L'Hopital's 
rule in Eqs. H27|l and (|28|l which after some lengthy calcu- 
lations gives 



A 



dp- 



1 2 



dR 



B 



d0^ 
dR 



C = 



A = 1 + (r - 1 



dlnPl 



B 



d\n. p 
d\n(5l 



d\n p 



(4 - ip) 
*8(1-m) 



C = 



(4-3/^)3 



1 - 



d\n p 



(F-l) 



(4-3/^)2 



(1-m) 



(31) 

(32) 
(33) 
(34) 



As we are interested for wind type solutions, the 
derivative dfP/dR\^ at the sonic point should be posi- 
tive. By noting that A,B > 0, this is equivalent to the 
condition C < 0. Then, from Eq. (|34|l and employing eqs. 
1)181) . H19() . and H29I) we get an inequality for /t. We solved 
this numerically and found the condition p > pmin where 
the value of pmin as a function of F is shown in Fig. [3 as 
the lower limit of the grey filled region. 

Thus, the criticality conditions for wind- type solutions 
reduces the acceptable values of /z to a limited range: 



Pn 



< p< Pn 



(35) 



Mmax = '"G/''*,min (sce Fig.|21) corrcsponds to the mini- 
mum distance between the Mach surface and the black 
hole, r*^niin, for which Eq. (|18|l has an acceptable solu- 
tion p* = p*iPs*)- For the hmiting case r* — > r^,min, 
Pi, — > CO. If the sonic surface is too close to the black 
hole, gravity becomes so important that the critical 
point disappears. 



^ Pmin — ''G/''*,max (sce Fig.|21) correspouds to the max- 
imum distance, r* max, between the sonic surface and 
the black hole above which the two critical solutions 
are accretion-type solutions. 

From Fig. |21 we see that there is a maximum value, 
F = 3, beyond which /imin > Mmax, because the cooling 
is too strong to have a wind type solution similarly to 
Begelman (1978). 



0.6 - 



0.4 



0.2 



y =3 



sill 

Jii^ ^7=1.01 



7 =1.04 





1.5 



2.5 



Fig. 2. Contour plots of constant terminal Lorentz factor 
7oo are shown within the limiting values of /i, Pm\n and 
/imax versus F, with the grey filled region showing the al- 
lowed wind space parameter. Note that /tmax corresponds 
to 7oo — > oo. 



3.3. Bernoulli constant 

The conservation of the number of particle can be inte- 
grated to give a first constant of the system, which corre- 
sponds to the mass loss rate or the mass flux m 



a'jnvr 



m . 



(36) 



Similarly, by integrating Euler's equation Q, we ob- 
tain the Bernoulli equation E — ajw, which is given by 



1/2 



1 



i + (i + r)p/2 



(37) 



Note that asymptotically (i? — > oo, p — > 0) the latter 
equation implies that the Lorentz factor is 700 = E/cq- 
On the other hand, by using the values of /3* = Ps* and 
P-i, at the sonic surface 

(eqs. j2ni and ^HDj we get the Bernoulli constant 700 — 
E/eq as a function of F and p, see Fig. 13 

Instead of integrating the differential equations which 
is the way we followed after Parker (1960), we could also 
plot the isocontours of constant E and m. The two meth- 
ods are equivalent. We have checked that while integrating 
the differential equations E remained constant. 
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4. Results and comparison with the nonrelativistic 

cases 

In the following we present the results of our parametric 
study for various values of the polytropic index T and the 
gravitational parameter fi. We also compare them with 
the corresponding results in the analysis of Blumenthal 
& Mathews (1976) for an adiabatic polytropic index T = 
5/3 and studies using a classical equation of state with a 
constant polytropic index T. 

4.1. Parametric study 

We solved numerically the equations by using a Runge- 
Kutta scheme starting at the sonic surface 

and integrating both upstream towards the black hole 
and downstream to infinity. 

The behavior of the solutions for various values of fi is 
displayed in Figs.|2| As /i increases the sonic surface 

moves inwards closer to the compact object. Thus, the 
supersonic region of the outflow exposed to thermal en- 
ergy deposition is enlarged. Consequently the Bernoulli 
energy gets larger for a given gravitational energy. This 
induces a better conversion of thermal energy into kinetic 
energy. In other words, as /i increases the distance be- 
tween the corona and the black hole decreases and, to 
support gravity, the pressure becomes higher. The pres- 
sure driven wind passes through the critical point faster 
and subsequently transforms into a relativistic wind at 
large distances from the black hole. The situation is 
similar to that of nonrelativistic winds with momentum 
and heat addition (,Leer fc Holzer 1980) wherein only for 
heat/momentum addition in the supersonic part of the 
flow the terminal wind speed increases. The mass loss is 
unaffected in such a case and is only affected for momen- 
tum/heat addition in the subsonic region. As illustrated 
by Fig.l^b, the higher is fi, the more extended downstream 
is the area of the flow where the effective polytropic index 
is smaller than F = 1.2 and thus the wind is more effi- 
ciently accelerated. This more efficient conversion of ther- 
mal energy for higher values of fi maybe also illustrated 
by calculating the cross-section of the equivalent De Laval 
Nozzle UParker 1960|l . 



din 5, 



'eff 



dR 



2 
R 



1 - /jg 



(38) 



where ScS is the effective cross-section. As shown in Fig. 
13;, the larger is /i, the larger is the opening of the effective 
cross-section. 



4.2. Comparison with the classical wind 

In the classical limit, the horizon goes to the center of the 
gravitational well, re —^ and consequently /i — > 0. The 
index F is the only free parameter which is left together 
with the mass loss rate. 



Then, in the classical limit the condition at the critical 
point reduces to the usual expression. 



PL 1 2 SM. 
^ 4 2r* 



(39) 



where the sound speed at the sonic surface 

is equal to half of the escape speed. From this single 
relation, the distance between the stellar surface and the 
sonic surface is uniquely determined from the sound speed 
at the sonic surface 

Conversely in the relativistic case ^ is a free parameter, 
so we can change the distance between the two surfaces 
and increase the effect of thermal conversion in accelerat- 
ing the wind. 

For a classical polytropic wind there is no acceleration 
for F > 3/2 fsee lParker 1960|) because then the accelera- 
tion at the sonic surface 

becomes negative for both critical solutions and no real 
wind-type solution exists. Conversely our relativistic solu- 
tions can be accelerated even if F > 3/2 when /i increases. 
This is illustrated in Figs. 0Ji,b where the topology of the 
adiabatic solution (F — 5/3) is displayed in the classical 
and the relativistic cases. In the first case both critical so- 
lutions are decelerated while in the second one there is an 
accelerated wind type solution. This is also shown in Fig. 
121 for /i = (classical limit, along the horizontal axis), the 
highest value of F to have an accelerated wind-type solu- 
tion is 1.5. For non zero values of (relativistic regime), 
the domain of wind-type solutions extends to higher val- 
ues of F. 



The relativistic wind is characterized by the fact that from 
the compact object outwards the pressure decreases faster. 
This is due to the presence of the redshift function a in the 
continuity equation (|36(l . When matter flows away from 
the central object the proper volume of the fluid increases 
faster than the observed (spherical) volume because of a, 
which provokes the decrease of the proper density and 
consequently of the pressure. The larger pressure gradi- 
ent leads to higher acceleration. This is equivalent to the 
well known acceleration due to overradial flaring in the 



classical solar wind as studied in Kopp fc Holzer (1976) 

In fact, expanding Euler's equation l|^V|l to first order 
in the relativistic effects, i.e. for <C 1 and p <C 1, we find 



1 (/3//3s 

2 1 - 



1 91n/32 



2i?2 (32 



/32 

F(F 



dR 



i?2 



(40) 



The first term of the right hand side is the classical term. 
The second term shows the effect of gravity in general rela- 
tivity, related to the existence of a characteristic distance 
due to the metric of the black hole. This positive term 
increases the acceleration and the conversion of thermal 
energy into kinetic energy. The third term shows the influ- 
ence of relativistic temperatures. It is also a positive term 
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Fig. 3. For F = 1.2 we show respectively as functions of the radial distance, the Lorentz factor in (a), the effective 
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Fig. 4. Comparison of the classical Parker solutions and the relativistic ones for two sets of the parameters F and fi. 
For an adiabatic polytropic index F = 5/3 and /i = 0.65 in a) is shown the nonaccelerating classical Parker solutions 
while in b) the accelerating relativistic ones. For F = 1.2 and /i = 0.2 is shown in c) the classical Parker solutions and 
in d) the relativistic ones. 



which increases the efficiency of the acceleration as long as 
the temperature remains relativistic. These two terms are 
responsible for the acceleration obtained even for F > 3/2. 

Even for quasi-isothermal (F 1) flows we see that 
the relativistic winds are more rapidly accelerated than 
their classical counter part (Fig.^,d). 

4.3. Comparison with a relativistic wind having a 
constant effective polytropic index 

Let us consider next relativistic winds with a constant ef- 
fective polytropic. In this approach the same free param- 
eters exist, /Lt and F. However for a relativistic polytropic 
wind where Feff = F =const, the expression of the sound 
speed and the Bernoulli equation are written as 

/3j = r-£-^, p..^, (41) 

(F - 1) + Fpc neo 



E = a-fi^l + ^r-^Pcj eo . (42) 

The main difference in the behavior of the solutions with 
those using the generalized equation of state appears in 
the zone close to the compact object where there is a tran- 
sition between classical and relativistic temperatures. 

For a fixed value of fi, we may compare the solutions 
of three equations of state describing a relativistic wind. 
First, by using the consistent equation of state with a 
variable Fcs and second, for a constant polytropic index, 
either Fj-ei = (F-|- 1)/2, a value which is the limit of Feff in 
the ultrarelativistic domain, or, for a constant Fnomoi = F 
which is the limit of Foff in the classical regime (see Fig.lSJ. 
The consistent relativistic equation of state gives a wind 
which is always more efficient than the classical one cor- 
responding to the constant value F, but less efficient than 
the ultrarelativistic one for the value (F -I- 1)/2. Using the 
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value r will always underestimate the asymptotic speed, 
but this difference becomes smaller as we approach the 
adiabatic value, as illustrated in Fig. |31 Conversely, it is 
appropriate to use the value (F + l)/2 in the vicinity of 
the black hole where the temperatures are ultrarelativis- 
tic. However it will always overestimate the asymptotic 
Lorentz factor, especially in the adiabatic case. 

As it can be seen in Fig. |SJ3, in some cases the solu- 
tions can be slightly decelerated in the super-sonic region, 
because the thermal energy is not sufficient to overcome 
gravity in these distances. 



the relativistic polytrope is more efficiently accelerating 
the flow than the classical one if F > f . 

At this point of the comparison, we should stress that 
a constant polytropic index, as it has been used in many 
models of relativistic thermal wind, is inconsistent with 
having both ultrarelativistic and classical temperatures in 
the flow. With the usual equation of state the thermody- 
namic regime is either ultrarelativistic or classical. In re- 
ality the outflow escapes from a very hot corona but cools 
down further out so there must be a smooth transition 
from one regime to the other. 




Fig. 5. Comparison between a variable and two constant 
F polytropic equations of state for two solutions with /i = 
0.3. In the left panel F — 1.2 while in the right panel F = 
1.6. In each panel, the two constant values Frci = (F-|-l)/2 
and Fnonrei = T corrcspond to the limits of Fcs in the 
ultrarelativistic and the classical regimes, respectively. 

In the ultrarelativistic limit (large fi and p) we can 
obtain analytically that our model is more efficient than 
the classical one, by comparing the acceleration at the 
sonic point for the two models. The condition is 



dR 



< 



dR 



(43) 



where the subscript C refers to the classical polytrope, 
and R to the relativistic one. This is fulfilled if 



1-/3; 



s \c 



1 



s \C 



21- fi 
which is equivalent to 



> 2 



1 



(3^1 1 
\r ^ 



(3f 



s \R 



2 1 



(44) 



(45) 



This condition is satisfied if the sonic surface 

is close enough to the compact object. In the limit of 
large ^, comparing Eqs. (|1T|| and (fll^ . we conclude that 



5. Application 

In this section we apply our model using typical values 
for jets from compacts objects. We consider two cases, 
a supermassive black hole with M = 10^ M©, typical 
of a quasar (see |Wang fc Zhang 2003| ) and a stellar size 
black hole with M = 10 Mq, typical of a microquasar (see 
IMirahel IflflSj i. 

We choose the sonic surface to be close to the compact 
object, Ti, = 6.45 T-Q, such that ^ = 0.155. This choice 
of /i gives a sufficiently high asymptotic Lorentz factor 
of 7 > 5, typical of those objects for sufficiently small 
values of F (Fig. EJ. For the solution with F = 1.09, the 
initial temperature is T « lO^^K in the case of electron- 
proton gas and T k, lO^K in the case of electron-positron 
pairs. Note however, that for low values of F such as F = 
1.09 the temperature remains unrealistically high at large 
distances. Hence, a more physical approach would be to 
avoid taking these solutions corresponding to F ^ 1 all 
the way to large distances, but instead match them with 
solutions with a faster decrease of the temperature (F « 
4/3). 

In both cases the solutions differ only in the density 
and mass loss rate. Both solutions have the same Lorentz 
factor profile and the same energy per particle in units 
of mass energy as this depends only on /i and not on the 
mass of the central object or on the total mass loss rate. 
The result is displayed on Fig. |B| 

For quasi-isothermal winds, F 1 — 1.2, the Lorentz 
factor 7 can exceed a value of 5 if the thermal energy dis- 
tribution in the vicinity of the central object is roughly 
10 times more than the mass energy. The thermal con- 
tent is spread along the flow and not peaked close to 
the compact object such that it can accelerate efficiently 
the flow with a lower maximum in temperature. Such a 
quasi-isothermal outflow can occur if there is a redistribu- 
tion of energy by the highly radiative initial field which 
cools the upstream flow and reheat the down part in the 
subsonic part where the medium is still optically thick 
(Das 2000; ITarafdar 19881 . The dissipation of disorga- 
nized magnetic field could also occur to produce extended 
heating. Such a dissipation has shown to be potentially 
important after the sonic surface and in the asymptotic 
region ( |Heinz fc Begelman 2000| l 

In Fig. [3 we plot the number density in units of 
no cm~^ as a function of distance, for a free density pa- 
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Fig. 6. Plot of the Lorentz factor in (a) and the thermal 
energy density in (b) versus radius for /z = 0.155 and three 
different values of the polytropic parameter F. 
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Fig. 7. Plot of density versus radius for /i = 0.155 and 
F — 1.09 for the quasar solution in (a) and the micro- 
quasar solution in (b). 



rameter uq- For the quasar solution (Fig. [7^) the corre- 
sponding mass loss rate is Mquasar = 10"^no MEddington = 
1.41 X lO~^noM0/yr, while for the micro-quasar solu- 
tion (Fig. |7|d) M;.quasar = 10" ^"^rio MEddington = 1-41 X 

10-22^0 Mo/yr. 

We note that as the mass of the central compact ob- 
ject changes from lO^Af0 to lOAf© the solution is simply 
scaled down spatially, a result consistent with the idea 
that microquasars could be considered to zeroth order as 



scaled down versions of quasars. A more detailed treat- 
ment however, should take into account certain differences 
of the two cases, e.g. pair production, different densities, 
etc., something beyond the scope of the present analysis. 

6. Conclusion 

We have generalized a variable polytropic index equation 
of state for the purpose of modelling relativistic flows, 
both in temperature and velocity in the vicinity of a 
Schwarzschild black hole. This has enabled us to analyze 
thermally driven winds having both ultrarelativistic tem- 
peratures at the base of the central corona (fceT <; eo) 
and classical temperatures (A;bT <C eq) further out. This 
equation of state is characterized by a polytropic index F 
such that pressure is related to density in the form of Eq. 

For a given polytropic index F, transonic wind solu- 
tions can be found only within a limited range of radii of 
the sonic point, , 

1 1 
< — < . 

This sonic transition should be far enough from the 
Schwarzschild radius tq such that gravity is not too high 
to allow the existence of critical solutions, because for too 
strong gravity the flow cannot escape, even if the temper- 
ature is infinite. On the other hand, the sonic transition 
should also be close enough such that the corona is not too 
diluted and the pressure too low to obtain an accelerated 
transonic solution. 

Schwarzschild's metric tends to enhance the effects of 
gravity. One major effect of strong gravity is to have a 
more efficient De Laval nozzle which allows to have accel- 
erated winds even for polytropic indices larger than the 
typical Parker's value, i.e. F > 3/2, conversely to the clas- 
sical Parker's wind IjParker 1960|l . This can be understood 
by computing the effective polytropic index c? In P/c? In p. 
Enhanced gravity and also relativistic temperatures tend 
to lower the effective polytropic index in the low corona 
which gives a more efficient thermal driving of the wind. 

However, we note that despite its widespread use in 
the literature, the ordinary polytropic equation of state 
with a constant effective polytropic index seems not to be 
really consistent with a mixed regime of temperatures in 
the corona from ultra- relativistic to classical ones. 

In order to reach very high Lorentz factors, in the adia- 
batic case, the internal energy of the plasma in the corona 
should exceed by a large factor its mass energy. This can 
be easily achieved if the total pressure is not limited to the 
kinetic pressure of the gas but also includes extra physical 
processes such as MHD waves or radiation. This result is 
somehow consistent with the usual low velocities obtained 
for classical winds. Quasi-isothermal winds reach higher 
Lorentz factors with a relatively lower - but still larger 
than the mass energy - internal energy. 

Applying our model to typical values observed in mi- 
croquasars and quasars we recover as expected that micro- 
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quasar outflows may be seen as a scaled down version of 
their bigger extragalactic counterparts. We do not claim 
however that all jets are only thermally driven. A magnetic 
driving mechanism seems more efficient indeed to acceler- 
ate the outflow on parsec scales as shown by Vlahakis & 
Konigl (2004), in agreement with some jet observations 
IjSudou et al. 2000)l . However, we simply emphasize that 
thermal driving may indeed play an important role along- 
side other mechanisms, when hot coronae are observed; 
in such cases, a consistent way of dealing with these rela- 
tivistic temperatures is required. 

Finally, this consistent generalization of the Parker 
polytrope for relativistic thermal winds could be imple- 
mented in numerical simulations, instead of the classical 
constant polytropic index equation of state to simulate the 
transition from relativistic to nonrelativistic temperatures 
along the flow. Nevertheless, as in any polytropic equation 
of state the source of heating is not specified on physical 
grounds and more detailed physics of the coronal heating 
is needed. 
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Appendix A: Heating of polytropic flows 

Conversely to the usual polytropic equation of state used 
in studying stellar interiors, the use in classical as well as 
relativistic winds of a polytropic index different from the 
adiabatic one is an implicit way to mask extra heating in 
the corona as explained in [Sauty et al. (1999) In all cases, 
the fluid remains a monoatomic plasma of electrons and 
protons with a ratio of specific heats 5/3. 

Thus, by substituting Eq. lO to Eq. (O, we have 
introduced implicitly that the external medium could give 
extra energy to the fluid, e in Eq. (|14|l is the total internal 
energy of the gas. It has two component, i.e. epiasma the 
internal energy of the plasma itself calculated from Eq. 

m 



which is in our case given by 



-plasma - ^0 ^ = ^5 + '° ^ > 

and the heating integrated along the streamline from the 
source Q. From a different point of view Q is also the 
internal energy (defined within an arbitrary constant of 
course) of the external medium (also called in thermo- 
dynamics "universe" as opposed to the system) which is 
given to the system through heating. 

The global system plasma + external medium 
being isolated, it is adiabatic, so it is usual to 
define e instead of Cpiasma in order to keep its 
simple form to the energy conservation given 
by Eq. Q (e.g., IBlumenthal fc Mathews 19761 
[Mobarry fc Lovelace 1986| ). It is easy to calculate 
the heating transfered to the plasma during its expansion. 



(A.3) 



The same holds for the enthalpy, w is the total en- 
thalpy of the gas and the external medium and if Wpiasma 
is the real enthalpy of the gas itself. 



W = Wpla 
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